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The idea of controlling the dynamics of spatially extended systems using a small number of localized 
perturbations is very appealing - such a setup is easy to implement in practice. However, when 
the distance between controllers generating the perturbations becomes large, control fails due to 
increasing sensitivity of the system to noise and nonlinearities. We show that this failure is due to 
the fact that the evolution operator for the controlled system becomes increasingly nonnormal as 
the distance between controllers grows. This nonnormality is the result of control and can arise even 
for systems whose evolution operator is normal in the absence of control. 



PACS numbers: 02.30.Yy, 05.45. Gg 

Control of spatially extended systems has recently 
emerged as a problem of fundamental interest as well as 
significant technological importance. Numerous investi- 
gations have shown the possibility of controlling a system 
by applying feedback at every point in space. However, 
the practical implementation of such control algorithms 
is often prohibitively expensive and sometimes impossi- 
ble. As a result, increasing attention is being devoted to 
localized control, where feedback is applied only at a few 
spatial locations. Previous studies have shown that the 
minimal number of such control locations depends both 
on the symmetry of the system |jj and the strength of 
noise [||. 

The failure of localized control for large separation be- 
tween controllers was attributed to the phenomenon of 
nonnormality in the evolution operator For nonnor- 
mal systems, stability becomes a poor predictor of short 
term dynamics p|. Strong nonnormality, which makes 
the system extremely sensitive to noise, was previously 
identified as the mechanism that provokes the transition 
to turbulence in uncontrolled systems, such as pipe or 
channel flows m ■ The latter systems are nonnormal due 
to a large mean flow, whereas generic extended systems 
have no mean flow and are normal for typical bound- 
ary conditions. In these normal systems, nonnormality 
arises as a result of control. Studies of nonnormality 
caused by localized control || have so far been limited to 
the interplay between nonlinearity and the noise ampli- 
fication due to nonnormality, rather than the emergence 
of nonnormality itself. Our goals here are to investigate 
how nonnormality arises in an originally normal system 
and study how nonnormality leads to noise amplification, 
which determines the limits of localized control. 

Let us consider the Ginzburg-Landau equation (GLE) 



(x,t) = <j>{x,t) + 4>"{x,t) - 4> 3 {x,t) 



(1) 



as a prototypical system. Although simple enough to 
allow analytical treatment, the GLE describes a very 
generic reaction-diffusion system. We therefore expect 



the main results of the following analysis to apply to 
most extended dynamical systems of this type. 

The unbounded system (|l|) possesses an unstable uni- 
form steady state, = 0. Our control objective is to 
make this state stable. The symmetry of an unbounded 
system requires at least two independent controllers in or- 
der to control a uniform target state [Q , so we will break 
the symmetry by imposing certain boundary conditions. 
Specifically, we will require that ip vanishes on one of 
the boundaries, e.g., 4>(0,t) = 0, reducing the minimal 
number of controllers to one. The only controller will 
be placed at the opposite boundary, x = I. With this 
arrangement the length I of the system plays the role of 
the distance between multiple controllers in a system of 
larger size. We choose the feedback law to be of the form 



<f>'(l,t) 



K(y)0(y,t)dy, 



(2) 



such that the control signal depends on the state (f> of 
the system on the whole domain (this condition can be 
easily relaxed). The feedback gain K(y) describes how 
each point inside the domain contributes to the feedback. 

The spectrum of the unperturbed linearized system is 
discrete with eigenvalues and eigenfunctions given by 



A„ = 1 - q%, f n (x) = sin(q n x), 



(3) 



where q n — (n — l/2)ir/l and n = 1, 2, ■ • •. These eigen- 
functions are normal, as expected, and form a convenient 
basis for the stability analysis of the perturbed system. 
Dropping the cubic term in ([!]) and projecting the re- 
mainder onto the basis {/«} we obtain 

oo 

K = AA - (-1)" Km ® m = ( M$ )"> ( 4 ) 

m— 1 

where $ n and K n are the Fourier coefficients of (f> and K 



4>{x) sui(q n x)dx, 



2 f l 

K n = - I K(x) sm(q n x)dx, 
' Jo 



(5) 
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FIG. 1. Strength of feedback as a function of the size / of 
the system. The maximum of if (x) in real space is compared 
with the Fourier coefficients ifi and K2 calculated using M . 
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The control problem for the PDE (Q) is thus reduced to 
finding an infinite set of coefficients K m , which will make 
the matrix M stable. 

It turns out that the structure of M simplifies the prob- 
lem remarkably. Suppose we take annxn truncation of 
M by discarding all rows and columns except the first n. 
Setting the eigenvalues of the truncated matrix, M n , to 
a sequence A^ , A 2 , • • • , A^ is then equivalent to solving a 
system of n equations linear in K m 's. In particular, one 
can change the first m — 1, • • • , s eigenvalues from A m to 
\' m and leave the rest unchanged by setting 



(7) 



for m < s and zero otherwise. Since the right-hand- 
side does not depend on the size of the truncated matrix 
M n , the result also holds for the full matrix M, so we 
can drop the index n. This is a very important result, 
because it allows us to calculate Fourier coefficients of any 
stabilizing feedback gain. It also allows us to determine 
how these coefficients scale with the size of the system 
I. Substituting (||) into (0) and after some algebra we 
obtain 



7T 



(2m-l)n; =1 (A m -A;) 



1)1(8 



(8) 



To make the matrix M stable we only need to change 
the first s positive eigenvalues, where s is equal to the 
integer part of l/ir+ 1/2. In particular, in the limit X[ = 
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FIG. 2. The magnitude of the disturbance in the middle 
of a large domain, I — 25, under LQR control. Exponen- 
tial asymptotic decay is preceded by a transient, amplifying 
the initial disturbance (white noise with standard deviation 
a — 10 -10 ) by six orders of magnitude. 



A' 2 = • • • = A' s = A, where A is some negative number, 
the product in the numerator of (||) reduces to (A m — 
A) s . Expressing s through I we see that the coefficients 
K m grow exponentially fast with I. The leading order 
behavior for large / is given by 



K„ 



exp 



I I / TYl 

- 2 + log(A m - A) + O -5- 



(9) 



This exponential growth is clearly seen in Fig. |. Both 
here and throughout the paper we use and compare two 
different control methods. In linear-quadratic (LQR) 
control the system (Q) is truncated to 64 Fourier 
modes and a set of K m that minimizes a quadratic form 
in $ is sought numerically. In pole placement (PP) the 
feedback gain is calculated directly from (j8|), where we 
change all unstable eigenvalues to A = —0.5 and leave 
the stable ones unchanged. (The choice of A is some- 
what arbitrary and is chosen to roughly correspond to 
the average of the first s eigenvalues produced by LQR) . 
Both control laws show the same scaling of K with I. 

The exponential growth of the control signal suggests 
transient behavior, a sign of nonnormality. Indeed, a 
small initial disturbance inside the domain will create a 
large control perturbation at the right boundary, x — I. 
If the feedback is designed properly, this perturbation 
will eventually (after propagating through the system) 
cancel the initial disturbance, thereby making the system 
asymptotically stable. However, asymptotic decay will be 
preceded by a transient whose magnitude grows with the 
feedback gain K . Numerical simulation of the linearized 
GLE does indeed show a large transient (see Fig. ||). 

The scaling of both the transient amplification and the 
control signal can also be understood qualitatively. Since 
the propagation of perturbations is diffusive, it will take 
a time r roughly proportional to the size of the system, 
I, for the control signal to travel from the right bound- 
ary to the left one, suppressing the disturbances inside 
the domain. Since the system is locally unstable, a dis- 
turbance near the left boundary will grow uncontrolled 



2 



during this time interval. The exponential growth of the 
disturbance will result in its amplification by a factor 
exp(Air) ~ exp(/3Z). To suppress the amplified distur- 
bance, we need to apply a control perturbation at least 
as large as the disturbance itself, which requires expo- 
nential (with I) growth in the feedback gain. 

Next, let us consider how nonnormality arises in our 
system. The evolution matrix for the controlled dynam- 
ics can be written in the form M = A + BK^ , where A 
is the diagonal (and hence normal) matrix that describes 
the dynamics of the unperturbed system, A mm = A m , 
and B and K are vectors with elements (— l) m + 1 and 
K m , respectively. Clearly BK^ is not a normal matrix, 
so the sum A + BK^ is not normal either. In fact, when 
all unstable eigenvalues are chosen equal, the matrix M 
is not even diagonalizable. In that case, one can convert 
M into the Jordan normal form 



J(M) = S^MS = 



Ji 
J 2 



(10) 



with Ji being an s x s Jordan block with eigenvalue A, 
J2 a diagonal matrix with A s +i, \ s +2i ■ • • on the diagonal, 
and S the respective coordinate transformation. This 
means that the eigenvectors of M corresponding to the 
Jordan block J\ coincide, e\ — &i = •■■ = e s . The 
solution of the linear system (||) in this case is constructed 
using the generalized eigenvectors e' p , such that Me' p = 
e' p _ 1 for p = 2, • • • , s and e[ — e\. Specifically, 



$(*) = 



E 



E 



Cp-\-m 



p—1 \_m— 



exp(Ai)e p 



^ c p exp(A p t)e p , 

p=s+l 



(11) 



where c±, C2, • • • are integration constants which have to 
be chosen to satisfy the initial condition. The result for 
the individual Fourier coefficients can be written more 
conveniently using the elements of the transformation 
matrix S: 



p=l 



s—p 



ECp+m . 



m=0 



exp(Ai) 



^ S n , P c p exp(\ p t). 

p=s+l 



(12) 



By looking at ( pL l| ) or (12) one can clearly see that the 
solution 4>{x, t) grows as a polynomial of order s— 1 before 
exponential decay at a rate A finally takes over. 

In the above analysis we required the first s eigen- 
values to be equal. What happens in the typical case 
when these eigenvalues are close, but not equal? First 
of all, as (JsJ) shows, the coefficients K m of the feedback 
gain will still grow exponentially with the size of the sys- 
tem. The matrix M will remain strongly nonnormal, 




FIG. 3. First five eigenfunctions of the controlled system 
for Z = 25. The evolution matrix M was calculated using 
LQR and has all eigenvalues distinct. 



but will become diagonalizable. Therefore, all eigenvec- 
tors of M will be distinct, but the first s will be closely 
aligned. This can be argued in the following way. Since 
K m quickly grows with increasing |A^J, while the size 
of the largest possible perturbation is usually severely 
restricted by practical limitations as well as nonlincar- 
ities, the new eigenvalues have to be chosen in a small 
strip of negative values, say (— e, 0). Therefore, the in- 
crease in the system size will force progressively more 
eigenvalues to lie in this strip, making the difference be- 
tween successive eigenvalues shrink at least as fast as 
e/l. As we have seen previously, setting s eigenvalues 
equal produces an s x s Jordan block, which causes s 
eigenvectors to merge. Since s is arbitrary, such merging 
will occur for any number of eigenvectors corresponding 
to identical eigenvalues. As a result, pairs of successive 
eigenvalues will continuously approach each other, align- 
ing the respective eigenvectors. (The continuity can be 
checked by a straightforward application of perturbation 
theory.) Fig. ^ shows that already for I = 25 the first five 
eigenfunctions are very closely aligned: all have a nearly 
sinusoidal shape with either 8 or 9 nodes, i.e., almost co- 
incide with the first stable (unperturbed) eigenfunction 
fg(x) = sm(8.5irx/l) which has 8 nodes. 

Finally, it is useful to derive the quantitative result for 
transient amplification, 



7 = max ■ 
t,*(o) 



l*(*)lls 



max(|exp(Mt)||2, 



(13) 



because this is the ultimate measure that determines 
when modal analysis and linear control break down. Let 
us again assume that all unstable eigenvalues are made 
equal and analyze the structure of the solution ( |l2| ) more 
carefully. The transient occurs because each of the terms 
(t m /ml) exp(Ai) first grows as t m and then decays as 
exp(At), reaching the maximal value at t m = —m/A. 
This maximal value is given by 



{t m y 



■ exp(At r 



1 



(14) 



v / 27rm|A| m 

so for small |A| the term with m = s — 1 dominates fll^ ) 
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A good approximation for the transient amplification 7 
can, therefore, be obtained by picking the initial state 
<I>(0) that corresponds to setting c m = 5 m . s and calculat- 
ing the maximum of the ratio ||$(t)||2/||<I>(0)||2: 



7 



= 1 °n,l t 

T^T m t ax 



-l e Ai 



(«-!)!' 



(15) 



where we only keep terms with p = 1 . It should be noted 
that the chosen initial state corresponds to a "near opti- 
mal" disturbance, rather than the "optimal" disturbance 
producing the largest transient amplification j| . 

If the matrix S is normalized such that S±. s — 1, it can 
be shown that, for arbitrary n, S n s = 5 n 1 and 



c -I iy +n (^-A) 8 nL 2 (At-A) 
Since all S^i scale in the same way, ([Lq) gives 



7 



|A|« 



(16) 



(17) 



In order to perform the calculation of the leading order 
behavior of S^i we assume that A = 1 — Up — l/2)ir/l] 2 , 
where p is some integer, in which case ( |l6|) gives 

(A 1 -A)- 1 ( P + S -1)! 



\Si,i 



(P- 



l)\p{p- l)s!(s- 1)!' 



(18) 



We can now re-express p in terms of A thus extrapolating 
between the known expressions for A's corresponding to 
integer p's. At the leading order we again obtain expo- 
nential growth with the system size: 



7 ~ exp 



M l0gA A 



1 



log 



(a- l) c 



(19) 



where a = \\ — A. 

As Fig. |J shows, the numerically calculated transient 
amplification factor does indeed grow exponentially fast 
and is rather insensitive to the way the feedback gain is 
calculated. The slope is seen to be slightly different from 
the one predicted by ([l9|). This is to be expected. First, 
we only consider the terms in (^|) for which m = s — p. 
While that gives a correct leading order result for small 
|A|, for A = —0.5 the contribution from m = s — p — 1 is 
about half that from m — s — p and more terms might 
need to be considered. Second, ( |l9| ) is the asymptotic re- 
sult for ( jllf ) and is valid only in the limit of large I (here 
for I > 20). On the other hand, the numerical accuracy 
in calculating the matrix norm decreases rapidly with /. 
This is in fact a numerical fingerprint of nonnormality. 
The results found using standard numerical routines are 
getting rather inaccurate for strongly nonnormal matri- 
ces, e.g., for large / (in our case also for / > 20). 

The large transient amplification makes the system ex- 
tremely susceptible to noise, as noise is amplified by feed- 
back before being suppressed. In order for linear control 



■ Vpp 

° \qr 
si -Eq.(17) 
■- Eq.(19) 
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FIG. 4. Transient amplification 7 as a function of I. The 
squares show the values calculated numerically using ([TJ|) . 



to work, the magnitude of the nonlinear terms has to be 
smaller than the magnitude of the linear terms. Com- 
parison of their relative magnitude can be used to esti- 
mate when noise will start to interfere with control. For 
instance, for cubic nonlinearity in the GLE, using the 
argument of Egolf and Socolar || we obtain that the 
largest magnitude of noise, a, tolerated by linear control 
should scale like 7 -3 / 2 . However, our numerical calcu- 
lations produce different scaling. This disagreement is 
currently under investigation. 

To conclude, we have shown that the application of 
spatially localized control to spatially extended systems 
renders the system nonnormal. The degree of nonnormal- 
ity quickly increases with the size of the system due to the 
close alignment of a progressively larger number of orig- 
inally unstable eigenfunctions. Increasing nonnormality 
leads to transient amplification which quickly grows with 
the size of the system, thus imposing strict limitations on 
the density of controllers required to control a system of 
a given size in the presence of noise or truncation errors. 
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